#!/bin/csh
set d=1
set n1=128
set n2=128
set beta=0.
set ki=4.0
set famp=0.1
set ftype=deterministic
set nu=0.001
set mu=0
set tmax=0.5
set scheme=rk3
set snapshots=4
set courant=0.1
set inicond=KH
set Np=0
set Nsp=0
#-----------------------
switch ("$ftype") 
    case deterministic : # deterministic forcing 
        echo 'deterministic forcing'
        breaksw 
    case stochastic    : # stochastic forcing 
        echo 'stochastic forcing' 
        set famp=-$famp
        breaksw 
    default:
        echo 'no forcing ?'
endsw
#-------------------------
echo " -----------------------------"
echo  "Device =$d"
echo "grid=$n1 X $n2" 
echo "QGS beta = $beta,viscosity = $nu, Ekmann term = $mu"
echo "Forcing, amplitude = $famp, at knumber =$ki"
echo "tmax=$tmax"
echo  "time stepping scheme $scheme ..."
echo "... with courant number = $courant" 
echo "store $snapshots snapshots"
echo "Initial condition is $inicond"
echo " -----------------------------"
./ihd -d $d $inicond -f \ $famp $ki -$scheme \
   -b $beta -m $mu -n $nu  -t $courant $tmax -s $snapshots $n1 $n2 \
   -P $Np -S $Nsp 
